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The response of spherical two-spin interaction models, the spherical ferromagnet (s-FM) 
and the spherical Sherrington-Kirkpatrick (s-SK) model, is calculated for the protocol of the 
so-called nonresonant hole burning experiment (NHB) for temperatures below the respective 
critical temperatures. It is shown that it is possible to select dynamic features in the out-of- 
equilibrium dynamics of both models, one of the hallmarks of dynamic heterogeneities. The 
behavior of the s-SK model and the s-FM model in three dimensions is very similar, showing 
dynamic heterogeneities in the long time behavior, i.e. in the aging regime. The appearence 
of dynamic heterogeneities in the s-SK model explicitly demonstrates that these are not 
necessarily related to spatial heterogeneities. For the s-FM it is shown that the nature 
of the dynamic heterogeneities changes as a function of dimensionality. With incresing 
dimension the frequency selectivity of the NHB diminishes and the dynamics in the mean- 
field limit of the s-FM model becomes homogeneous. 



PACS Numbers: 64.70 Pf,05.40.+j,61.20.Lc 



I. Introduction 

Non-exponential relaxation behavior is found to be rather common when dealing with dis- 
ordered materials like glasses, spin-glasses, disordered crystals or proteins pQ. In the last 
decade particular attention has been payed to the question to which extent the relaxation 
is to be viewed as dynamic heterogeneous [2]. Different experimental techniques have been 
invented in order to investigate the detailed nature of the relaxation particularly of amor- 
phous systems |S1 in El El- These techniques allow to specifically select a slow sub-ensemble 
and afterwards monitor its relaxation. Various applications of these techniques have shown 
that the primary response in amorphous polymers and supercooled liquids is to be viewed 
as heterogeneous in the sense that it is possible to select slow sub-ensembles relaxing at 
smaller rates than the average. Throughout this paper I use the definition given in ref.jZj 
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according to which a system will be called dynamic heterogeneous if it is possible to select 
dynamically distinguishable (slow and/or fast) contributions to a relaxation. 

In addition it has been found that after a certain re-equilibration time the relaxation 
properties of the selected sub-ensemble return to those of the bulk|8n9j. Therefore, these 
experiments indicate that the response can be described as a superposition of exponentially 
decaying entities with different relaxation rates. The various relaxation rates, however, 
are not static quantities but apparently fluctuate in time. Different interpretations have 
been provided for this behavior jTHl HI] • The NMR-techniques|Hlin] have the advantage of a 
simple interpretation in terms of equilibrium 4-time correlation functions but are restricted 
to a rather narrow temperature regime and to certain materials. The optical deep bleach 
technique j4| has the advantage that it can be applied in a wider temperature range with the 
shortcoming that up to now it has not been interpreted in terms of equilibrium correlation 
functions. At this point it is important to note that these experimental methods monitor 
molecular reorientations and therefore are not able to address the question of spatial aspects 
of the dynamic heterogeneities. The only exception is provided by a new variant of the 
4d-NMR technique jl2j, which allows to extract a length scale via monitoring spin-diffusion. 
Originally, this technique has been applied to a polymeric liquid and later on the length 
scale of the dynamic heterogeneities has been extracted for low-molecular glass-forming 
svstems fT^ ITl] . The corresponding length scales have been found to be on the order of 
1 ■ ■ ■ 4nm. 

Another experimental technique allowing to monitor dynamic heterogeneities is pro- 
vided by the nonresonant spectral hole burning (NHB) experiment |15j. This method is 
based on a pump-wait-probe field sequence with a large pump field amplitude beyond the 
linear response regime. Though originally applied to supercooled liquids ^Hl; NHB in the 
meantime has been used to investigate the relaxation of several materials, including dis- 
ordered crystals (relaxor ferroelectrics) ^/Tj , ion-conducting glassesjjlSj and spin glasses [T^. 
Also the application of NHB to a solvable glass-model has been presented |2()j. The interpre- 
tation of the obtained results has mainly been guided by the fact that via the application 
of a large amplitude ac field of frequency Q the sample absorbs energy of an amount 
proportional to the imaginary part of the susceptibility evaluated at the pump frequency 
f2|2l]. In case of a homogeneously broadened response one does not expect that it is pos- 
sible to modify the response in a frequency selective way. By contrast such a goal could 
be achieved if the response is given by a superposition of differently fast relaxing entities 
(heterogeneous scenario). This is because in this case energy absorbtion should be largest 
for those sub-ensembles with a relaxation time on the scale of the inverse pump frequency. 
This intuitive picture is confirmed in the framework of a response-theory for NHB for the 
particular case of stochastic dipole reorientations, which I have developed recently p2|. 

In the quoted experiments it always has been found that a frequency selective modifica- 
tion of the response indeed is possible. Regarding the re-equilibration of this modification, 
however, the results differ not only with respect to the frequency dependence but also 
regarding the time scale of the recovery. For the latter a time scale longer than the inverse 
burn frequency has been observed in the case of the relaxor ferroelectrics|17j . 
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As long as one is concerned with supercooled liquids, one can safely consider the system 
to be in (metastable) equilibrium prior to the NHB field sequence. This, however, is not 
necessarily true for the relaxor materials or the spin glasses. In particular, it was argued 
in ref.j2n| that the results obtained for the solvable p-spin-glass modelj2ni are mainly to 
be interpreted in terms of out-of-equilibrium effects. In equilibrium, the p = 3-model stud- 
ied shows an exponential relaxation at long times. Therefore, according to what was said 
above one does not expect to be able to select a sub-ensemble in a frequency dependent 
way in the equilibrated version of the model. It has to be mentioned here that dynamic 
heterogeneities as monitored by NHB have also been observed in Monte Carlo simulations 
on an equilibrated Sherrington-Kirkpatrick mean- field spin-glass model|2Sj. These calcu- 
lations explicitly demonstrate that from the observation of dynamic heterogeneities one 
cannot conclude on the existence of spatial heterogeneities. Also, this finding appears to 
be independent of whether the system was in thermal equilibrium before the application 
of the pump-field. 

In this paper I consider the application of the NHB field sequence to the spherical 
Sherrington-Kirkpatrick (s-SK) model, i.e. the p-spin model with p = 2 and the spherical 
ferromagnet (s-FM) in arbitrary dimension. The Langevin dynamics for these models has 
been solved analytically [23 • Additionally, the s-FM model is equivalent to the 0(N)-model 
in the limit of large N and therefore is a typical model for domain coarsening processes |26j. 

I will solely consider a thermal history protocol in which the system is quenched to 
a temperature below the critical temperature from infinite temperature prior to the 
experiment. The behavior at and above will be investigated in a forthcoming publication. 
Therefore, all effects observed are intimately related to the out-of-equilibrium dynamics 
of the model. This is because the system never reaches equilibrium in this temperature 
regime. The outline of the paper is as follows. In the next Section I will briefly recall the 
dynamic features of the models and calculate the response to the NHB fleld sequence in 
second order regarding the pump-fleld amplitude and linearly in the small step-fleld. In 
Section HI the results of the calculations are presented and discussed. The paper closes 
with some conclusions in Section IV. 



II. NHB in spherical two-spin interaction models 

The spherical models under consideration are deflned by the Hamiltonian 

H = JifcSjSfc - ^/ijSj (1) 

where in case of the s-SK model the Jik are chosen at random from a Gaussian probability 
distribution with zero mean and variance a = 1/N and are restricted to a ferromagnetic 
coupling J (to be set to unity in the following) for the s-FM model on a simple hypercubic 
lattice in d dimensions. In addition the spin-variables are subject to the spherical constraint 
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Y,i s~ = N. The Langevin equations governing the dynamics of the model read as 

Siit) = E JikSkit) + hi- z{t)si{t) + ii{t) (2) 

k 

where z{t) is the Lagrange multipher enforcing the spherical constraint and ^j(t) is a 5- 
correlated Gaussian white noise. These equations have been solved analytically by Cuglian- 
dolo and Dean (CD)|23] in two papers. The dynamical properties of the s-FM model are 
discussed in ref. j^Zl- Further information regarding the correspondence between the two 
models can be found in ref . [2H] • Here, I briefly summarize the results relevant in the present 
context. 

Of particular importance are the violations of the fluctuation-dissipation theorem 
(FDT), which relates the response function to the time-derivative of the two-time cor- 
relation function, 

In particular, it has proven extremely useful in out-of -equilibrium situations to define the 
so-called fluctuation-dissipation ratio X{t^tw) via|29j. 

T3(+ + \ X{t,tu,) dC{t,t^) 
R{t,t^) = —^ 9^ 

the limiting long-time behavior, X^, of which is known to vanish for domain coarsening 
modelsjSO]. Here, denotes the time that has elapsed after a quench to the working 
temperature prior to the measurement. Also for the models considered in the present 
paper one has = 0^7\. Thus, concerning this measure of typical distances from 
equilibrium, the domain coarsening models in finite dimension do not show any differences 
to the mean-field s-SK model. 

The response of the system, Rh{t,t') = Yli^ {si{t)^i{t')) /{2NT), in the presesence of a 
field h{t) can be obtained in the same way as calculated by CD for the zero field case. As 
shown by Berthier et al.fHT] this yields: 

.,-0^.(4^) (5, 

where the function g{t) is defined by 

g{t) = [exp (-4t)/o(4t)]'' s-FM 

^(t)=exp(-4t)^ s-SK (6) 
with In{x) denoting the generalized Bessel function and Wh{t) is the solution of 

Whity = g{t) + 2T fdTWh{rfg{t - t) 
Jo 

+ J^dh J^dt2h{U)h{t2)Wh{ti)WH{t2)g(t - ^-^^'^ (7) 

which has its origin in the normalization of the equal-time correlation, C{t,t) = 1, i.e. the 
spherical constraint. 
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Zero field response 

Before I turn to the calculation of the response following the NHB pulse sequence, it is 
appropriate to summarize the known results for the reponse and the correlation in zero 
field, for a more detailed discussion see refs. [23 123 I2H1- In zero field, eq.© simplifies to 
the following Volterra equation: 

W{tf = g{t) + 2T fdrW{rfg{t - r) (8) 
Jo 

For the s-SK model, this equation has been solved by CD for random initial conditions 
(i.e. a quench from T = oo at t = 0). No simple analytical solution exists for the s-FM 
model. However, for t ^ 1 it can be shown [2^1 EZl ESj that for T < Tc 

Wit)' = -^^—l^g,,{t) (9) 

Here, the asymptotic behavior of the g{t) are given by 

gasit) = (87rt)-^/2 s-FM 

= (327r) -1/2^3/2 s-SK (10) 

The critical temperatures are given by Tc = 1 for the s-SK model and Tc depends on the 
spatial dimension in case of the s-FM model, with Tc{d = 3) ~ 3.9568jn2]- For the other 
dimensions used in the present paper I find Tc{d = 5) ~ 8.6482, Tc{d = 7) ^ 12.7982 and 
Tcid = 9)- 16.8579. 

If tw denotes the time elapsed after a quench from T = oo to T < Tc, the following 
behavior is found for R{t + t^, tw) from eqns. ()5p9pi()|l in the hmit of long t^: 

R{t + tw, tw) = X-'/^g{T/2) with A = ; > 1 

tw + T 

e = d (s-FM) ; e = 3 (s-SK) (11) 

From this expression it is evident already that the dynamic properties of the s-SK model 
are very similar to those of the s-FM in c? = 3. Therefore, the following expressions are 
given for the s-FM model. The only differences between the s-FM model in c? = 3 and the 
s-SK model stem from the different prefactors in eq. (fTn|) . In particular, it is evident from 
eq. ljlll) that the temperature is an irrelevant variable in the whole low-temperature phase. 
Two time-sectors are to be distinguished: 

• For short times r such that t <^tw - the so-called stationary regime - one has A ~ 1 
and accordingly: 

R{T + tw,tw)=R{T)=g{T/2) T<^tw (12) 

In particular, in this stationary regime the FDT, eq.(jHI), holds. In the long time limit, 
1 <^ T, eq. lfT^ shows that R{t) decays according to R{t) ~ r"'^/^. 
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• In the so-called aging regime 1 ^ r ~ t^, one finds from eq. (jll|) : 



"W 



(13) 



In this regime the FDT is strongly violated. If in addition r ^ t^, the response 
behaves as i?(r, t^,) ~ r"'^/'^. 

The overall behavior is shown for various values of the waiting time in Fig.l. From 
this plot it is seen that for all dimensions shown, d = 3,5, 7, there is a crossover from the 
r"*^/^ to the T~^^'^ behavior. Also the explicit dependence on the waiting time in the aging 
regime is evident. It should be pointed out that the crossover from the stationary regime 
to the aging (domain growth) regime takes place around r ~ independent of spatial 
dimension. Remember, that only in the stationary regime the FDT holds. In particular, 
the behavior of the response (and also the two-time correlation function) does not change 
qualitatively around d = 4, above which the model behaves mean-field like concerning the 
statics and the exponents. (The same holds for the dynamic fluctuations, as will be shown 
elsewhere.) 

The above discussion shows that the relaxation is extremely non-exponential. There- 
fore, the question as to what extent the response can be viewed as dynamic heterogeneous 
naturally arises. 

Nonresonant hole burning 

In the following, the response will be calculated for the NHB-field-sequence, cf. Fig.2. At 
t = the system is quenched from T = cx) to the working temperature T. One (or more) 
cycles of the pump-field hp{t) = hpsin [Q{t — tg)] are applied after a time tg has elapsed. 
Following a waiting time the response R*{i + T,i + r') is measured, where I defined 



for brevity. Here, N denotes the number of cycles of the sinusoidal pump-field. In the 
following calculations of R* the Wh{t) are needed in second order with respect to the 
pump-field amplitude hp. From eq.(|7j) it is evident that a perturbation expansion follows 
from 



Inserting this expression into eq.(|7)) yields eq.® for the zero'th order term and, assuming 
a time-dependent field of the form h(t) = hp sin [Q{t — tg)]9(t — tg) according to Fig.2: 



t = tg + tp + and tp = 2N7C / Q 



W,^{tf = W{tf + hlA{tf + 0{h'^) 



(14) 




tl+t2 

2 



) 



(15) 
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where =Min(r, tp). Using eq. (|14p the response function in 0{hp) then is found to be 
given by (r > r'): 



R* {i + r, £ + t') = R{i + r, £ + r') + Ai?(£ + r, £ + r') 

^ Vr(£ + r') (t-t'\ 

R{t + r,t + T')= , ^^ (71 * 



Vr(t + r) 



(16) 



Ai?(£ + r,£ + r') = 



A(£ + r)2 A(£ + r 



/\2 



l^(t + r)2 vr(t + r')2 



i?(£ + r, £ + r') 



In the NHB-protocol of Fig.2, however, the response to a small step field is recorded, i.e. 
the integrated response (the thermoremanent magnetization): 



X*(£,r) = x(t,r) + Ax(£,r) = / dsR*{i+T,i+s 



(17) 



according to eq. lfTBj) with the zero-field integrated response x(^)T) = JodsR{i + r, £+ s). 
Eqns. fll5|) - ■ -(^Zj) allow the calculation of the results of a NHB experiment at any desired 
temperature. 



III. Results and discussion 

Many of the general features of the modification Ax(£, r) can already be seen for T = 0, 
which is the simplest case. Afterwards finite temperatures will be discussed as well as the 
dependence of the observed features on spatial dimension. Finally, a direct comparison 
between the three-dimensional s-FM model and the s-SK model will by carried out. 

Throughout the remainder of the paper the dependence on £ will be skipped whenever 
no confusion can occur, i.e. the shorthand notation x(r) = x(^, t) and Ax(r) = Ax(£, r) 
will be used. Times and frequencies will be given in dimensionless units. 



A. Spherical ferromagnet 

In this subsection the details of the results for calculations of the response following the 
NHB field sequence are discussed. The discussion is kept general with regard to spatial 
dimension but the actual calculations are carried out for d = 3, cf. Figs. 3- ■ -6. The depen- 
dence of the results on spatial dimension will be presented in the next subsection. 



T=0: 

For T = 0, eq. (fTH|) can be solved trivially and the modified response is easily calculated. 
From eq.(jHI) one explicitly has 

W{ty = git) and /\{tf = Ao{t) ; T = (18) 

The corresponding expressions for x(t) and Ax(r) are easily obtained from eqns. ()16|17|) . 
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From the discussion of the zero field response in the last section it is evident that the 
time tq elapsed after the quench and before application of the pump-field is a very important 
parameter. For small tq some transient features are expected due to the interplay of the 
approach of the aging-regime and the additional non-equilibrium features induced by the 
application of the pump-field. Of course, the cross-over to a f2 independent behavior will 
depend crucially on the pump-frequency Q as this determines the time tp of the imposed 
non-equilibrium situation. This is demonstrated in Fig. 3a, where Ax(t) is plotted for 
d = 3, T = 0, = 0, Q = 0.1, 10 and several values of tq. From this figure two features 
are evident immediately. First of all it is seen that Ax(r) is non-zero only in a limited 
time interval and the time of the maximum modification depends on the burn frequency 
fl, thus demonstrating dynamic heterogeneous behavior. Additionally, the curves for tq = 
differ from the others in that they change from positive to negative values in a limited time 
range. This transient behavior also depends on Q. However, it is always possible to choose 
tq in a way that the mentioned interplay between the two sources of transient features can 
be neglected. Thus, it is interesting to consider the asymptotic regime determined by 

tq 1 where tq = tq{Q). 

The detailed dependence of tq on Q has to be found empirically in the sense that no transient 
effects should show up in the modified response for a given pump frequency Q. In order 
to further demonstrate the relative independence of the results from the chosen value of 
tq, in Fig.Sb I have plotted the value of the maximum modification, Axmax = ^xi'i'max) 
vs. tq for Q = 0.01,0.1,10 in a scaled way. This plot demonstrates the features already 
mentioned above. For small tq there is some time interval in which Ax(r) is negative 
and therefore the maximum is reduced. The small hump around tg ~ 27r/f2 in the curves 
is roughly located at those values of tq where also the negative part vanishes. Finally 
for long tq a plateau is observed and the results are independent of tq. Obviously, the 
minimum value of tq satisfying the constraint that Ax(r) is independent of tq depends on 
Q. In the whole range of Q considered in the present paper it turned out that a value of 
tq = 10^ is sufficient. Therefore all further calculations are performed for tq = 10^ unless 
stated otherwise. (I have checked via explicit calculations that xitq + tpj ^) ~ xi^q: r) ~ 
independent of f2. Otherwise, the interpretation of the f2-dependence of Ax(£, r) in terms 
of dynamic heterogeneities would be meaningless.) From eqns. ljlOp and ()15|) it is seen that 
in this regime one has for Ao(tg + r) for arbitrary dimension [tm =Min(r, tp)): 

Ao(t, + r) = (Sn)-'/' fj^ dt, fj'dt^ sin {Qt,) sin {nt^) ^^^^ ^^[^^^^^ ^ ^^^j,^, » 1 (19) 

This expression along with eq. (jl6|) for the modification of the response also explains the 
observed relatively weak dependence of Ax(r) on tq. In lowest order Ao(t + r') behaves 
as t""^^^ and according to eq.Q the same holds for W{i + r'). Therefore, R and Ai? are 
independent of tq in this order and the same holds for A^. An explicit t^-dependence enters 
only in higher order. 
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For T = 0, the above expression is only needed for tm = tp = 2Ntt/Q, where denotes 
the number of cycles of the sinusoidal field. The reason for the fact that smaller values of 
are irrelevant in this case is that according to eq. lfT^ ^(t)^ = ^o(^) and that the response 
is measured only after the pump-field is switched off. Therefore, the question as to which 
extent the results depend on the number of cycles naturally occurs. Fig. 3c shows Ax(r) as 
a function of r for = 1, 5, 10 for Q = 0.1. The curves for = 5 and = 10 are hardly 
distinguishable in this plot. The saturation of the maximum amplitude is demonstrated in 
Fig. 3d, where Axmax is plotted vs. A^ for = 0.1 and 10 (tq = 10^). For A^ values larger 
than roughly 6, Axmax become A^-independent. In all of these calculations one has tg S> tp, 
whereas the opposite limit is met at small tq as discussed above, although for A^=l. This 
demonstrates that the field sequence of NHB is unable to drive the system much further 
away from equilibrium than it already is due to the quench at t = 0. 

After this consideration of the influence of the parameters tg and A^ now the more 
important issue of the fi-dependence of Ax(r) will be discussed. Before presenting the 
results of model calculations it is instructive to investigate analytically the limiting behavior 
of Ao(tq + r) according to eq. ()19|) which in turn determines the behavior of Ax(r). One 
finds that Aq vanishes in the limits of large and small Q according to 

Q^O: Ao(t + r) ~ ^^^^ ^ ^-^ 

Ao(t + r) ~ {Qtg-^l) (20) 

n^oo: Ao(t + r)~^]-^ Vd 

demonstrating that the modification induced by the NHB pulse sequence in principle de- 
pends on Q. The universal fi~^-dependence for large Q can easily be understood by slightly 
rewriting eq. ()15|l in the form {M = Qtm) 

A„((, + r) = £'dt, £'dt, s,u s,u t,W (t, + |) («, + |) j (. - ^i-t^i) 

The behavior for large Q is obtained via second-order expansions of W{tq + Q~^x) and of 
g{T — Q~^x). In a simple calculation one then finds the quoted behavior. Therefore, 
this represents a universal result which is valid for any model obeying eq. ()15j) . As the Q- 
dependence of Ao(tg + t) uniquely determines the one of A{tq + tY, this result additionally 
holds for all temperatures, including the disordered paramagnetic phase. 

In Fig. 4a the dependence of Ax(r) on Q is demonstrated in detail for d = 3. It is clearly 
seen that Ax(r) shows a very pronounced f2-dependence for Q < 0.1 which, however, 
diminishes with increasing Q. Additionally, it is evident that the spectral modifications 
become 'broader', i.e. are non-vanishing in a larger time interval, as Q decreases. 

In Fig. 4b I have plotted the time of the maximum modification, Tmax vs. Q. From this 
plot one can see that Tmax varies as for small Q. This is the behavior typical for an 
extremely broad distribution of relaxation times. Additionally, r^ax becomes independent 
of Q for Q ~ 10. Thus, for the s-FM model in d = 3, the response is dynamic homogeneous 
in the short time regime, whereas it is dynamic heterogeneous for long times. Therefore, 
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the aging dynamics in this model is dynamic heterogeneous. Also, it is seen from Fig.4b 
that the behavior is independent of the time tq elapsed after the quench. 

Another important question regards the time scale of the 'recovery' of the modification, 
i.e. the waiting time dependence. In ref. [221 I have found that in a model of reorient at ional 
dynamics there is no extra time scale for the recovery. However, in the experiments on 
the relaxor materials a very long recovery time scale has been found^Tj. In order to 
investigate this question, in Fig. 5a I have plotted the maximum modification (normalized 
to the value at ty^ = 0), Ax{Tjnax)norm., vs. scaled waiting time fltu, for a variety of burn 
frequencies Q. It is evident that the life-time strongly depends on Q for small Q and that 
this dependence diminishes for larger Q. It has to be mentioned at this point that the form 
of the modifications hardly change for different waiting times. In order to have a simple 
measure for the recovery times, I have fitted the curves of Fig. 5a to a Kohlrausch function, 
exp { — (ttf,/rR-)^^}, and plotted the resulting time scales tk vs. Q in Fig. 5b. I included 
the time of the maximum modification for = (cf. Fig. 4b) in this plot (dot-dashed 
line). As the behavior of these two time scales is almost identical to within a factor of two, 
the conclusion to be drawn from these calculations is that there does not exist an extra 
time scale for the recovery in the s-FM model. This appears plausible on intuitive grounds 
because the only time scale set by the pump is roughly 1 / {27cQ) in the domain growth 
regime and just unity (l/(27r)) in the short-time regime. This is qualitatively different 
from a complex domain-structured system like the relaxor materials, where pinning effects 
play a dominant role. 

T>0: 

So far, I have considered T = solely. As already mentioned in the last section, the response 
in the ferromagnetic phase is independent of temperature. For the NHB field sequence, 
however, eqns. ()15|) and (fTBj) show that AR and correspondingly also Ax(r) do depend on 
temperature. The physical reason for this dependence is quite clear. Though the thermal 
noise is irrelevant for the linear response, this changes under NHB conditions due to the 
aligning effect of the pump-field. Here, thermal fluctuations tend to counterbalance the 
induced alignment. 

In order to investigate the importance of this effect, eq. ()15|l is solved numerically for 
finite temperatures for tg^l, i.e. using eq.© in the expression for Ao(t). (To ensure that 
the physical solution of the equation is met, I first performed calculations for very small 
T, in which case eq. (ll5|l can be solved in terms of an expansion in T/Tc.) The results 
of such temperature dependent calculations are shown in Fig.6a, where Ax(r) is plotted 
for various Q and for temperatures up to 0.5Tc using tg = 10^ and ty^ = 0. Two features 
are evident by inspection of that plot. First, the position of the maximum modification 
hardly changes as a function of temperature. It should be mentioned that also the shape 
of the Ax(r) does not change as a function of temperature. This means that no effect of 
'motional narrowing' is observable. 

The most prominent feature is the increasing intensity of Ax(r). In order to demon- 
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strate this behavior in more detail, the intensity at the time of the maximum modification, 
Ax{Tmax), is shown as a function of reduced temperature l—T/T^ in Fig. 6b for fi = 0.01, 100. 
This plot demonstrates a scaling behavior of these quantities. The exponents, however, 
are different for various frequencies and do not seem to have an obvious explanation. It 
would be interesting to further analyze the f2-dependence of Ax{Tmax), which, however, is 
beyond the scope of the present study. 

The conclusion to be drawn from these calculations is that the effect of temperature is 
seen mainly in the amplitude of the modification, but narrowing effects are absent. This 
means that the nature of the dynamic heterogeneities are not affected by temperature- 
effects. Therefore, in the following, I will concentrate on T = for simplicity. 

B. Spherical-ferromagnet: NHB for varying dimension 

In contrast to the s-SK model the s-FM model offers the opportunity to study the behavior 
of the observed dynamic heterogeneities as a function of spatial dimension. As already 
mentioned in Sect. II the fluctuation-dissipation ratio X^q vanishes for the s-FM model 
independent of spatial dimension. Additionally, it is known [27j that the usual static critical 
exponents take on their mean- field values for d>4 while some dynamical exponents depend 
on dimension for arbitrary d. The only hint for a change in the dynamic properties of 
coarsening models for d > 4 stems from the large N-model, where it has been shown 
that the aging contribution to the integrated response changes qualitatively for d > 4^^]. 
The two-time quantities in the spherical models considered here, however, do not show 
any signature of a change in behavior around d = 4. Therefore, in this subsection, the 
dependence of the behavior of Ax(r) and thus of the dynamic heterogeneities on spatial 
dimension will be investigated for T = 0. 

Fig. 7a shows a plot of Ax(r) versus r for d = 5 (upper panel) and d = 7 (lower panel), 
which is to be compared to Fig. 4a. (I do not consider even spatial dimensions like d = 6, 
because the calculations are much more involved for technical reasons.) In both cases the 
il-dependence is much weaker than in d = 3. Whereas there is some weak il-dependence 
for d = 5, such a dependence is hardly visible for d = 7. To quantify this diminishing 
dependence, in Fig. 7b I have plotted the time, Tmax, of the maximum modification as 
a function of fl. Whereas one has a dependence in the long-time regime for d = 3, 
this is weaker, roughly fi"^/^, for d = 5 and finally in ci = 7 and d = 9 there is no visible 
f2-dependence. This means that the dynamics in (i = 7 an (i = 9 is dynamic homogeneous in 
the short-time and the long-time regime. Thus, also the aging dynamics becomes dynamic 
homogeneous for higher spatial dimension. 

The conclusion from these calculations is that the dynamics of the s-FM model becomes 
dynamic homogeneous in the mean-field limit, (i^l. However, as pointed out earlier, the 
mean-field limit holds for all (i>4 when concerned with static properties. For d = 5, Fig.7 
reveals that the aging dynamics still is heterogeneous, although the shape of an effective 
distribution of relaxation times is changed relative to d = 3. 
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C. NHB in the spherical Sherrington-Kirkpatrick model 

As already noted in the preceeding section the dynamic behavior of the s-FM model in 
d = 3 and the s-SK model is very similar in the low temperature phase. Concerning the 
modified response x*{'^) the same holds true, in particular in the limit of large tq. The 
only difference stems from the functions g{t) occuring in the expressions for Ao(tg + r) 
and R{i + T,i + r'), eqns. ()lfjpi9|l . for the two models. Plotting the two functions gpuit) 
and gsxit) reveals that gpAiit) — gsKit/'i) in a good approximation for moderate t and 
asymptotically one has gsKit) ~ [{Mf / {?>2T:Y/^]gFM{t) 12. hlgpuit), cf. For 
large Vt small times r are relevant and therefore one expects /\xsk{t) ~ Axfa/(t) to hold 
with the difference that the maximum modification occurs at — '^^max- ^^e other 
hand, for small Vt long times are most important and all functions g{t) can be replaced 
by their asymptotic values. Thus, one roughly has /S.xsk{j) ~ [(87r)^/(327r)] Axi?M(T) and 
^max — '^max- Fig. 8a I have plotted Ax(r) for both models for Q = 10"^ (upper panel) 
and Q = 100 (lower panel). Apart from a factor of two in the amplitude for the larger 
frequency the behavior just discussed is recovered. 

Fig. 8b shows the time of the maximum modification, Tmax versus frequency Q. The full 
line is for the 3d s-FM, cf. Fig. 4b, and the dashed line is for the s-SK model. Also included 
is (t^^/3) for > 1 (dot-dashed line). From the above discussion it is clear that one just 
finds the expected behavior: for very small Q one has ~ whereas ~ Sr^al 
holds for large fl. 

These considerations show that the qualitative behavior of the s-FM model for d = 3 
and the s-SK model are extremely similar in the low temperature phase, including the 
dynamic heterogeneities. The similarity between the 3d s-FM and the s-SK model is well 
known and it has been argued earlier that the s-SK model does not generically behave like 
a spin glas, concerning both, the statics and the dynamics j2H]- 

The fact that there are dynamic heterogeneities observable in the s-SK model might 
appear somewhat astonishing in view of the fact that the s-SK model is a mean-field 
model. However, this result is not unique. The first observation of dynamic heterogeneous 
relaxation in a mean- field model was reported in ref.j^Ol- Also the Monte Carlo simulations 
on the Sherrington-Kirkpatrick model mentioned above |21] show dynamic heterogeneous 
behavior. The present results therefore confirm the fact that an identification of dynamic 
and spatial heterogeneities is not possible in general. 

IV. Conclusions 

In the present paper I presented calculations of the response of simple spherical spin models 
to the field sequence of the nonresonant hole burning experiment, a technique allowing the 
detection of dynamic heterogeneities if they exist. The calculations were restricted to the 
generic situation of aging in the low temperature phase. It was assumed that the system 
is quenched from infinite temperature to T <Tc at the beginning, correlations in the initial 
conditions thus being neglected completely. The equilibrium dynamics will be considered 
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in a forthcoming publication. 

The main result of the present paper is the fact that the non-equilibrium dynamics 
of the s-FM model in ci = 3 is dynamic heterogeneous in the aging regime whereas it 
appears homogeneous in the short-time regime. In higher dimensions, d>7, the response 
becomes dynamic homogeneous also in the aging regime. Therefore, the nature of the 
aging dynamics changes and in the mean- field limit, d^l, the model displays homogeneous 
dynamics. This is what one would expect as in this limit each spin interacts with an infinite 
number of neighboring spins. However, for the s-SK model, a simple mean-field spin glass 
model, a behavior very similar to that of the s-FM model in d = 3 is found. This explicitly 
demonstrates that the existence of dynamic heterogenities does not tell us anything about 
spatial heterogeneities in the system. Furthermore, a strongly heterogeneous dynamics has 
also been observed in the aging behavior of a short range spin glass model [^1] showing 
strong relations between local responses and correlations. The mere existence of dynamic 
heterogeneities in a short range spin glass is of course to be expected and here also the 
spatial aspects of these heterogeneities are of importance. 

When considering the spherical two-spin interaction models, it is tempting to associate 
some disorder with the spherical constraint which forces the equal time correlation function 
to unity at all times. The lengths of the spins are not static quantities and are random 
to some extent. This might be viewed as a kind of 'dynamic' disorder and the similarities 
between the s-FM model in d = 3 and the s-SK model hint towards the irrelevance of the 
quenched disorder in the latter when concerned with dynamic quantities. On a speculative 
level the fact that the dynamics in the s-FM model becomes homogeneous for d ^ 1 
can be understood from the following argument. Assuming the existence of an effective 
distribution of relaxation rates at a given instant of time, i.e. for a given distribution of 
Si{t)^ the width of this distribution is expected to decrease with an increasing number of 
nearest neighbor interactions and consequently the dynamics becomes more homogeneous. 
This of course does not mean that the response decays in an exponential way for large 
d because the considered distribution is not a static quantity. Notice that this argument 
also implies that the life-time of the heterogeneities is finite. Unfortunately, NHB does not 
allow to determine this life-time. It is, however, important to point out that the argument 
concerning the increasing number of neighbors, though appealing with respect to the s-FM 
model, has important drawbacks. Taking it serious in case of the s-SK model would predict 
homogeneous dynamics in contrast to what is observed. 

Although the s-FM model is a typical model for phase-ordering kinetics, there is no 
obvious relationship between the observed dynamic heterogeneities and the domain size 
distribution in a coarsening system. A possible way to investigate such a relationship could 
be to perform calculations on the Ising model using spatially varying magnetic fields. It 
would also be interesting to perform an analysis along the lines of refs. j34| l35j in order to 
see whether the dynamic heterogeneities in a coarsening system behave similar to those 
observed in glassy systems. Furthermore, such calculations would allow to investigate the 
dependence on spatial dimension and therefore to check whether the one observed in the 
present paper for the s-FM model also is found for other domain-coarsening models. 
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In summary, I have shown that heterogeneous aging can be observed in the low- 
temperature phase of the spherical model of a ferromagnet. The aging dynamics becomes 
homogeneous on increase of the spatial dimension despite the fact that no qualitative 
change in the two-time quantities like the correlation function or the response is observed. 
Quenched disorder does not play any significant role with respect to heterogeneities in 
spherical models. 
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Figure captions 



Fig.l : The response function R{t + tw,tu,) for the s-FM model in the ferromagnetic 
phase, T < Tc, for tyj — 10^, 10^, 10^, 10^ (from upper to fowest hne). Upper panel: 
d = 3, middle panel: d — h, lower panel: d — 7. 

Fig. 2 : The field sequence for the nonresonant hole burning (NHB) experiment: A time 
tq after a quench from T — oo one or more cycles of a strong sinusoidal field 
hp{t) — hpSm{ft{t — tq)) are applied. After a waiting time the response to an 
infinitesimally small field is monitored. 

Fig.3 : a: The modification of the response, Ax(r) vs. r for the s-FM, d = 3. Upper 
panel: Q = 0.1; lower panel: Q — 10. The time after the quench from infinite 
temperature, tq, is chosen as tq — 0, 10, 10^. 

b: AXmax,sc. = ^X7nax(tq)/AXmax(tq^lO'^) VCrSUS tq for Q = 0.01,0.1, 10. 

c: Ax(t) versus r for f2 = 0.1 and tq — 10^ for iV = 1 (full hne), N = 5 (dashed 
line) and iV = 10 (dotted line). The curves for N — 5 and N — 10 are hardly to 
distinguish. 

= Axmax{N) / Axmax{N = 1) as a function of the number of cycles of 
the sinusoidal field, (tg = 10^). 

Fig.4 : a: Ax(r) vs. r for ^^ = 10^ t^^O, d^3 and log (O) = 1,0, -1, -2, -3, -4 (from 
left to right). 

b: Full line: Tmax vs. fl, dashed line: flTmax vs. fl for tq = 10^, t^ = 0. Additionally 
shown as the dot-dashed line is Tmax vs. fl for tq=l and t^, = 0, demonstrating that 
the behavior does not depende on tg. 

Fig.5 : a: Ax{T„iax)norm. vs. the scaled waiting time flt^j for fl = 10~^, 1, 5, 10, 25, 50, 10^ 
(from left to right). 

b: Characteristic decay time tk of a Kohlrausch-fit of the from Ax{Tmax)norm. — 
exp { — (t^/rx)'^^} {tg = 10^). The stretching parameter is approximately constant, 
— 0.9. The dot-dashed hne represents r^ax for tyj — for comparison, cf. Fig.4b. 

Fig.6 : a: Ax(r) vs. temperature for T = 0, 0.1, 0.2, 0.3, 0.4, 0.5 for the burn frequencies 
given in the respective panels. The remaining parameters are tg = 10^, tyj = and 
d=3. 

b: The value of the maximum modification, Ax{rmax), vs. 1 — T/Tc. The other 
parameters are the same as in a. 

Fig.7 : a: Ax(r) vs. r for burn frequencies Q = 1, 10, 20, 50, 100 for d = 5 (upper panel) 
and d—7 (lower panel). The remaining parameters are tg — 10^ and t^ — O- 
b: The time of the maximum modification, Tmax, for d — 3, d — 5 and d — 7 vs. fl. In 
the upper panel Tmax is shown and in the lower one the product flTmax- PuU lines are 
for tg — 10^ and the dot-dashed fines for tg — 0. 
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8 : a: Ax(r) vs. r for the 3-dimensional s-FM (full line) and the s-SK model (dashed 
line). Upper panel: Q = 10^'^ and the factors muliplying Ax(r) are N{SK) = 327r 
and N[FM) = (Stt)^. Lower panel: = 10^. The remaining parameters are tq = 10^ 
and tw — 0. 

b: The time of the maximum modification, r^as, for the 3d s-FM model (c.f. Fig.4b, 
full hne) and the s-SK model (dashed hue) vs. Q for the same parameters as in a. 
Also included is t^^/?> (dot-dashed line) for Q>1. 
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